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SUMMARY 


Woodward’s panel method for subsonic and supersonic flow is improved 
by employing control points determined by exactly matching two-dimensional 
pressure at a finite number of points. The results show great improvement 
in the predicted pressure distribution of a flapped airfoil. With the 
paneling scheme of cosine law in both chordwise and spanwise directions, 
the method is shown to accurately predict leading-edge and side-edge suction 
forces of various configurations in subsonic and supersonic flow. 

INTRODUCTION 

Woodward's unified subsonic and supersonic panel method (ref. l) with 
constant pressure panels has been in wide use for sometime. The predicted 
overall longitudinal aerodynamic characteristics by the method are reasonably 
accurate. However, the predicted lifting pressure distribution is not 
accurate enough for calculating the leading-edge and side-edge suction 
forces. The latter are quite sensitive to the paneling scheme and the 
control point locations. In reference 2, a paneling scheme derived from 
numerical experimentation was suggested for calculating the leading edge 
suction. Again, this approach is still not good enough for calculating the 
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side-edge suction force. Examination of the two-dimensional flat plate 
problem indicated that if a constant-percent control point location is 
used, the overall force and moment and the pressure distribution cannot be 
accurately predicted simultaneously. Originally, Woodward (ref. l) recommended 
95% of the panel chord as the suitable location for the control point. Later, 
85% location was suggested based on matching the two-dimensional lift of a 
flat plate (refs. 3,^)« However, 85 % location of control point does not 
produce accurate pressure distribution and pitching moment. 

Recently, Dillenius and Nielsen (ref. 5) used the panel method to 
calculate the leading-edge and side-edge suctions in supersonic flow. After 
the strengths of the panel singularity have been obtained, they replaced 
them by a system of equivalent horseshoe vortices. The in-plane forces 
are then calculated by applying Kutta-Joukowski theorem and assumed to be 
acting along the planform edges. Only few supersonic results have been 
compared with other known theoretical calculation. 

In this report, a new method to improve Woodward’s panel method will 
be described. The pressure prediction is improved based on two-dimensional 
theory. The method used in reference 6 is then applied to calculating the 
leading-edge and side-edge suction forces by directly using the predicted 
pressure distribution. The method is applicable in both subsonic and 
supersonic flow. 


SYMBOLS 


A aspect ratio 

C leading-edge suction parameter (see equation 23) 

C total lift coefficient 

Li 


C 


chord length 
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'ref 


"t 

h/ c 

■t -t + 
i ,J ,k 


v , Jte 


v ,se 


t 

x,y ,z 


cp 


z 0 (5) 


S 

c 


sectional induced drag coefficient 
sectional lift coefficient 

sectional pitching moment coefficient about local leading edge 
root chord 
reference chord 

sectional leading-edge suction coefficient 

sectional leading-edge thrust coefficient 

maximum camber height to chord ratio 

unit vectors along x-, y- and z-axes , respectively 

planform lift curve slope per radian at a = 0° 

leading-edge suction coefficient at one radian angle of attack 

side-edge suction coefficient at one radian angle of attack 

Mach number 


number of chordwise panels 
number of spanwise strips 

distance from the leading-edge along the camber line nondimen- 

sionalized with respect to the local chord 

side-edge suction force per unit length of tip chord 

wing-fixed rectangular coordinates with positive x-axis along axis 

of symmetry pointing downstream, positive y-axis pointing to 

right, and positive z-axis pointing upward 

x- coordinate of center of pressure 

x-coordinate of the leading-edge 

camber function 

angle of attack 
3z 

= tan (r~-) 

ot, 


b 


6 flap angle 

6 = 6 at the leading edge 

T accumulated vortex strength defined by Eq.. (28) 

y vortex density 

A leading-edge sweep angle 

X taper ratio 

£ - ( ) / c 

THEORETICAL DEVELOPMENT 

Airfoils 

Planar 

The main objective in the initial development is to predict exactly 
the pressure distribution of a flat-plate airfoil at those locations such 
that the lift and moment coefficients can be easily obtained by integrating 
the pressure. According to the Quasi-Vortex-Lattice Method (QVLM) (Ref. j) , 
these locations are best to be given by the cosine law distribution: 

L = (1 - cos ( (2k - l)rr/2N ) ) /2 , k = 1, ..., N (l) 

it. c c 

where N is the number of chordwise panels and chord length is assumed to 
c 

be unity. The chordwise paneling scheme will be based on the following 
cosine law: 


£ J+1 = (1 - cos(jTT/N c ))/2, j = 0, 1, ..., N c 


(2) 


Assuming the angle of attack to be one radian, the airfoil integral 


equation can be written as 
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1+7T = / 


AC ( £ ’ ) 1 
X 


K - V 


(3) 


If the chord is divided into N segments on each of which AC is 

c & P 

assumed constant, equation (3) can he integrated over each segment to give 
N 


Utt = 2 AC p U k Hn|(5j - ep/Uj+i - 5.)|, i = 1 N 


(k) 


where E,. is the control point location and E, . < E, < Equation (i+) 

will now be used to find E. such that the predicted AC = AC (£. ) 

i P k P k 

will be exact. Since the problem is a nonlinear one, it can be solved by 
the following iterative method. Differentiating equation (4) with respect 
to gives 


N 3AC 
_ c P, 


~ 3£. 

J=1 i 


£n 


Vlii 


N 


r ' _ f 1 ^ AC p ^ - £ 

^i+1 H j=l P k ^j+1 S 


° X (5) 


Solving equation (5) will result in a set of 3AC^/9^^ - values. If 

AC^ d ^ is the desired value and AC^'* is the computed one, an incremental 
P k P k 

change in control point location A£^ can be obtained from 


N 3AC 

2° AE. = AC (d) - AC (c) , k=l, 

i=l 3? i 1 p k p k 


, N, 


( 6 ) 


If the starting is chosen at 95% of panel chord, four or five iterations 

can produce essentially exact solution. Once AC are obtained, the lift 

P k 

coefficient is computed as 


c„ = f L AC dE = g- Z 77 AC sin6d6 “ \ 2 AC sin0. 

I op s 2op 2N ck=l P k k 


( 7 ) 
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Similar expression for can be derived. With the predicted pressure 

distribution, the leading edge suction can easily be calculated. 


Cambered 

The control points found by the above method can be used directly for 
cambered airfoils. In this case, the flow tangency condition to be 
satisfied is developed as follows. Let the camber shape be described by 


F(£,z) = z - z c (£) = 0 


( 8 ) 


The flow tangency condition is then given by 


dz . 

(V cosed + V sinoik + ui + wk) • (- -rrr i + k) = 0 

co co 


(9) 


or, in linearized form. 


» i r i &c p (s ' )dg ' 
V Utt o - 5 

(X) 


-sina + 


cosa 


dz 
c 


( 10 ) 


After AC is calculated by satisfying equation (10) , the aerodynamic 
P 

forces can be calculated by resolving in the appropriate directions the 

pressure force which is assumed to act normal and along the camber line. 

According to the incompressible Bernoulli equation, the pressure coefficient 

for lift at high a, C', is related to C at small a as follows: 

P P 


C’ 


P 


1 - [(V^cosa + u) 2 + (Vasina + w) 2 ] 
V 

CO 

( cosa) C 

1? 


-2( cosa)u 


( 11 ) 
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It follows that 


AC* = (cosa)AC 
P P 


( 12 ) 


where AC^ is calculated through equation (10). According to figure 1, 
the aerodynamic characteristics are given by the following equations: 


c n = f AC cosa cos(a - <5 )ds + c. sin(a — 6 ) 
l o p c t o 


= f AC cosa 
o p 


cos (a - 6 ) 

t — ~ d£ + c, sin(a - 6 ) 

coso t o 


(13) 


ca = / AC cosa 
i o p 


sin(a - 6 


cos<5 


d£ - c, cos (a - 5 ) 
t o 


(1U) 


c = - / £AC cosad£ 
m o p 


(15) 


-1 3z c 

where 6 = tan ( ) and s is the length of camber line nondimensionalized 

C oc, 

with respect to the chord, and 6^ is 6^ evaluated at the leading edge. 

The integrals in equation ( 13 ) — ( 15 ) can be evaluated by midpoint trapezoidal 
rule on the 0 - plane as illustrated in equation (T)* The above formulation 


can be easily shown to be consistent with the exact theory for a flat- 

3z c 

plate airfoil. For example, with = 0 in equation (10) the solution 

dt, 

- p 

obtained is the well known AC = U v(l-5)/£ sina and hence c = 27Tsin a. 

P t 

It follows that 


2 2 

c^ = 27Tsinacos a + 2wsin asina = 2Ttsina 
2 2 

c , = 27Tcosasin a - 27rsin acosa = 0 

d. 
i 

TT 

c = - — smacosa 
m 2 


(l6a) 

( 16 b) 


(l6c) 
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Flapped 

The method, described above can be extended to airfoils with flap 
deflection. In this case, the distribution of pressure points and panel sizes 
on the airfoil and the flap will be given separately by the cosine relation, 
following the QVLM (ref. 7)> For a flat airfoil with flap deflection 6^, 
the exact linear solution is given by (ref. 8); 


AC (£) = ha 
P 




Jin 


sin 


sin 



(17a) 


c^ = 2ira + 2<5 f (t + sinx) 


(17b) 


£.e. 


7T r /T . 

~ 2 01 ' f 5 + 


, sin( 2 t ) N 
SinT + r -) 


(17c) 


where 


cos (p =2^-1 


(18) 


COST = - 1 

and £ is the flap hinge location. Using equation (l7a) for the desired 
AC^ , the correct control points can again be obtained through iteration. 

It should be noted that in general the iterative method to find the correct 
control points for a flapped airfoil may not converge, as some control 
points may move out of the corresponding panels at some iterative steps. 
When this happens, that particular control point is reset at 95% of the 
panel chord before the iteration is resumed. The correct set of control 
points is taken to be the one with least square difference of ACp between 
exact and calculated distribution. 
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Wings 

For three-dimensional wings with straight leading edges, the chordwise 
distribution of control points and panels follows directly from the two- 
dimensional theory described above. The spanwise distribution of control 
stations and panel width is based on that used in the QVLM (ref. 7) in 
accordance with the cosine law. The control stations are given by 

y 1 = (b/2)[l - cos (in/ (N g + l))]/2, i=l , ..., W s ( 19 ) 

and the strip widths are given by 

y, - (b/2)[l - cos((2j - 1 )tt/ 2(.N + l))]/2, j=l, . N + 1 (20) 

J s s 

where is the number of spanwise strips. The inboard gap created by 

the scheme given by equation (20) should be eliminated, and the tip inset 
from equation (20) is retained, as described in reference 7- In general, 
the control stations will not pass through the panel centroids. From 
extensive calculation, the chosen spanwise paneling scheme appears to be 
the best for calculating the side-edge suction force. 

Once the pressure distribution is obtained, the interpolation by 
Fourier series is used to find the distribution of the leading-edge 
suction, the streamwise vortex density, y , and the side-edge suction. For 
example, to find the leading-edge suction parameter C, AC^ is first 
multiplied by — sin 9 and then developed in a cosine Fourier series: 

N 

T C 

f(9) = — AC sin0 = a + 2 a. cosj'0 
2 P o j 


( 21 ) 



i iii mi mi ■■ 1 1 i i in ■■■ 1 1 min ■■imi ii i 


where the Fourier coefficients a. can he computed by midpoint trapezoidal 

<3 

rule in terms of AC as follows : 

P 


a o = iF'o f(e)ae =h-2 f(e k> 

c k=l 


a. = — Z 71 f(0)cosj0d0 - 2 f(0 v )cosj0, 

J 7T O 


C k=l 




Near the leading edge, AC^ - C /(l-x)/x . Hence, 


lim “■ AC sin0 - C = a +2 a 


2 p 


0 j=i J 


The sectional leading edge suction is then given by (ref. 2) 


= (tt/8)C 2 (1-]^cos 2 A l e ) 1/ ' 2 /cos 2 A l e 


The side-edge suction per unit length of the tip chord is given by (see 
ref. 9) 


s (.x) = upG (x' 


where G(x) is defined by 


G( x) = /f- lim A - y/(D/2) 5 - y 

J V »/2 


To find Y , the conservation of vorticity is used: 
'x 
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_2L + _X 
3x dy 


0 


(27) 


If T(x,y) is defined as 


r(x,y) = - St Y (x',y)dx' (28) 

x £ y 

equation ( 27 ) can be shown to give 


= 3r(x,y) (29) 

Y x 3y ^ 

In the linear theory, y = AC /2. Hence, equation (28) can be written as 

y p 

r(x,y) = - £ ^~ f Q o | AC p sinGdQ (30) 

Using equation (21), equation (30) can be integrated in closed form. The 
differentiation in equation ( 29 ) is performed through the use of trigonometric 
interpolation formula. The detail can be found in reference 6. 

As can be seen from the above description, not only the leading-edge 
and side-edge suction forces can be predicted by the method, but also the 
distribution of y^ which is needed in calculating some lateral- 
directional stability derivatives. 

NUMERICAL RESULTS AND DISCUSSIONS 

Airfoils 

For a flat-plate airfoil in incompressible flow, all aerodynamic 
characteristics of interest can be exactly predicted by the present method. 
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The computed control points are in general located at Ql.5% for the leading- 

edge panel and moved to 97 .2% for the trailing-edge panel. For a cambered 

airfoil with or without flap deflection, only approximate results can be 

calculated. Figure 2 shows the comparison of the present predicted pressure 

distribution for a circular-arc airfoil with the results by conformal mapping 

(ref. 10). The lift coefficient and the center of pressure are presented in 

Table 1. It is seen that the present results are quite accurate. One 

indication of accuracy of the method is the magnitude of the calculated 

induced drag coefficient. The exact potential-flow solution is c _ =0. 

d. 

l 

For the circular-arc airfoil with h/c = 0.031U as shown in figure 2, the 

calculated c^ is about 0.1% of c^. It becomes 0.2% of c^ when h/c 
i 

is increased to 0.20. 

The predicted pressure distribution for a flapped airfoil with flap 

chord ratio of 0.30, 6 = 30° and a = 10° is presented in figure 3. 

Also shown is the result by the original Woodward's method, in which 

uniform panel size on the airfoil and flap is assumed. The calculated 

ACp is assumed to be at the panel midpoint. In both the original and the 

improved methods, seven and four panels are used on the airfoil and flap, 

respectively. The calculated c 0 and c are also presented in Table 1. 

m le 

The discrepancy in the calculated c 0 and c by the present method and 

X/ in 

the exact solution, as given in Table 1, is mainly due to the quadrature 

method used (i.e. midpoint trapezoidal rule) in integrating the pressure 

force. The method, which is also used in the QVLM, cannot account exactly 

for the logarithmic singularity in the AC distribution. Although the 

P 

original Woodward's method can predict c^ and c^ reasonably well (see 
Table l) , the calculated pressure distribution is inaccurate near the flap 
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hinge, as shown in figure 3. Note the present method and exact solution 
are in excellent agreement, to the extent, that the dashed curve cannot 
he distinguished from the solid curve. 

Wings 

In three-dimensional cases, it is assumed that the control point 
locations obtained from the two-dimensional theory are directly applicable. 
One important application of calculating leading-edge and side-edge 
suction forces is the prediction of vortex lift through Polhamus' method 
of suction analogy as extended in reference 9- According to this method, 
the total lift coefficient for a flat wing exhibiting edge-separated vortex 
flow can be written as 


2 2 
C T = K sinotcos a + (K „ + K ) sin acosa 

L p v,i6e v,se 


(31) 


where the lift factors K , K n and K are the lift curve slope 

p v,x,e v,se 

per radian and the leading-edge and the side-edge suction coefficients at 

one radian of angle of attack, respectively, from linear potential theory. 

To validate the present method, extensive comparisons will be made 

of the predicted K , K „ and K with other theoretical results. 

P v,xe v,se 

Furthermore, from convergence studies with the present method, it was 
determined that solutions converged rapidly with increasing number of panels. 
For example, with eight or nine chordwise panels and 10 to 15 spanwise 
strips , aerodynamic characteristics of simple planforms can be accurately 
calculated with confidence. Further increase in the number of strips 
will not change the calculated results significantly. The results to be 
presented below were mostly obtained with nine chordwise elements and 12 
spanwise strips, except those cases when the pressure distribution is desired. 


In these cases, 13 spanwise strips were used. Since the chordwise control- 


point locations are based on the two-dimensional incompressible flow theory, 

any deviation of the pressure distribution from the two-dimensional theory 

will lead to some inaccuracy of the method. This appears to be the case in 

predicting K for highly swept planforms in subsonic and supersonic flow 
v , He 

with leading-edge sweep angle >_ 60°. In both situations, the predicted 

K n was found to be slightly too high. This problem has been solved 
v , x-e 

empirically by moving downstream the control points of the leading-edge 
panels by a certain percent, A, of the elemental panel chord, where, 
for M < 1 , A = 0 for A <_ 60° 

A = 0.35 (A - 60) for A > 60° 

for M _> 1 , A = 3.0 for A <_ 60° 

A = 3.0 + 0.35 (A - 60) for A > 60° 


For cranked wings also, a uniform constant percent chordwise location for the 
first control point is used all across the span for simplicity. The sweep 
angle for determining its location is obtained by weighting the sweep angles 
with respect to the length of the leading-edges of each section. 

Typical pressure distributions in supersonic flow are presented in 

figures and 5 and are compared with the exact solution obtained from reference 

11. The most visible improvement in the predicted pressure distribution 

appears to be over the inboard section. To improve the predicted pressure 

distribution in the tip region, it was found that the number of spanwise 

strips has to be further increased. However, the inaccuracy in the 

predicted pressure distribution in the tip region apparently does not affect 

the accuracy of the overall force coefficients because of the decreasing 

chord. The lift factors K ,K 0 and K for the rectangular wing 

of aspect ratio 2.0 are presented in figure 6 as a function of Mach number 

and tabulated also in Table 2. The agreement of the present results with 
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Lamar's (ref. 9) in subsonic flow is excellent. This is particularly true 
in K and K 0 prediction in that the present results are indistinguishable 
from Lamar’s. The supersonic results also agree well with the exact linear 
theory (ref. 12 and 13). Note that the original Woodward’s method tends to 
predict higher K^. 

The results for a cropped delta wing of low aspect ratio are presented 

in figure 7 and Table 3» The agreement with other theoretical results 

appears to be reasonably good, except that the predicted ^ in subsonic 

flow are slightly higher than those given by Multhopp's method of Lamar 

(ref. 9)- (The effect of Mach number on K was first shown for these 

v,se 

two Wings in reference l4). 

The sectional leading-edge and wing-tip suction coefficients for 
cropped delta wing are plotted in figures 8 and 9 at Mach numbers 0 and 1.8. 

At both Mach numbers the sectional leading-edge suction coefficients agree 
very well with other methods (refs. 7» 9, 11 and 15) in the inboard region. 

The agreement is not as good in the tip region. With regard to the wing- 
tip suction coefficients, at subsonic speed the present method is in good 
agreement with that of Lamar (ref. 9)- At both subsonic and supersonic 
speeds the present method shows an overshoot near the leading-edge caused by 
the method trying to represent the exact continuous solution with finite 
panels. 

Additional comparisons for various configurations are presented in 
Table 4. The planform shapes used in comparison include straight, tapered 
wings and a family of double delta configurations of reference l6 which are 
illustrated in figure 10. The results indicate that the present method can 
predict edge suction force values that are in reasonable agreement with 


other methods , even for complex planforms . 
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So far only the lift factors and pressure distributions have been 
compared. The centers, of these forces- can also be accurately predicted. 
According to the method of suction analogy, the pitching moment of wings 
with edge -separated vortex flow can be calculated as (ref. 9)- 


C 


m 


r „ .2 Ze ^ „ .2 se-, 

LK smacosa — ^ + K „ sm a + K sm a J 

p c v,xe c v.se c 

14 ' j* 



( 32 ) 


The calculation of the center-of- force factors x , x„ and x is 

p &e se 

compared in Table 5 for the cropped delta wing of figure 7* The results 
indicate that the centers of edge suction forces are well predicted by the 
present method. 

CONCLUSIONS 


Based on the extensive comparison of present prediction with other 
theoretical results, it may be concluded that the present improved Woodward’s 
panel method is generally accurate in predicting the leading-edge and side- 
edge suction forces and the centers of these forces in subsonic and super- 
sonic flow. The good accuracy of the present method has also been 
demonstrated for cambered and flapped airfoils. Because of generality 
of the panel method, the present improved method can therefore be used 
not only to predict the vortex lift of complex planforms through the method 
of suction analogy, but also to calculate certain lateral -directional 
stability derivatives as well. 
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Table 1. Predicted Aerodynamic Characteristics for 
Cambered and Flapped Airfoils 


Airfoil 

Aerodynamic 

variable 

Exact 

Present 

Woodward 
(ref. 1 ) 

circular -arc , 

c p 

(ref. 10 ) 
0.941 

0.9345 

Hot Determined 

with h/c = 0.0314 
a = 0° 

X 

CP 

0.355 

0.35575 

M 

Flat Airfoil with 

c o 

(ref. 8) 
3.265 

3 .30 

3-212 

30 deg. Flap and 
Flap Chord Ratio = 

c 

-1.152 

-1-175 

-1.192 

0.30, a = 10° 

m Jle 

















Table 2. Comparison of Predicted Lift Factors for a Rectangular Wing of A 



K 

P 



! K 

v,36e 



K 

v.se 


Present 

Lamar 
(ref. 9) 

Woodward 
(ref. l) 

Exact 
(ref. 13) 

Present 

Lamar 
(ref. 9) 

Present 

Lamar 
(ref. 9) 

Exact 
(ref. 12) 

2.4783 

2.4763 

2.5744 

Wot Available 

1.4991 

1.4997 

1.5793 

1.5489 

Not 

Available 

2.4956 

2.4934 

2.5928 

ft 

1.5027 

1.5033 

1.5978 

1.5673 

It 

2.5507 

2.5479 

2.6514 

tl 

1.5136 

1.5142 

1.6582 

1.6272 

It 

2.6553 

2.6517 

2.7627 

Tt 

1.5313 

1.5323 

1.7804 

1.7483 

It 

2.8389 

2.8354 

2.9592 

If 

1.5537 

1.5558 

2.0292 

1-9939 

1! 

2.9811 

2.9783 

3.1126 

ft 

1.5634 

1.5665 

2.2791 

2.2354 

1! 

3.7943 

- 

4.0393 

If 

0. 

- 

2.5737 

- 

II 

3.6321 

- 

3.8114 

3-7575 

0. 

- 

1.8664 

- 

11 

2.9864 

- 

3.0975 

3. 0408 

0. 

- 

1.2647 

- 

If 

2.5304 

- 

2.6129 

2.5615 

0. 

- 

0.9929 

- 

1.0194 

2.2062 

- 

2.2725 

2.2262 

0. 

- 

0.8301 

- 

0.8507 

1.9630 

- 

2.0181 

1.9760 

0. 

- 

0.7190 

- 

0.7351 










Table 3- Comparison of Predicted Lift Factors for a Cropped Delta Wing 


A - 0.874, A - 0.4 and A = 63 ° 


M 


K 

. P 



K 0 
v,£e 



K 

v,se 



Present 

Lamar 
(ref. 9) 

Woodward 
(ref. 1) 

Present 

Lamar 
(ref. 9) 

Exact 
(ref. 12) 

Present 

Lamar 
(ref. 9) 

Exact 
(ref. 12) 

0. 

1.2968 

1.2789 

1.3632 

1.5297 

1.50U1 

Not 

Available 

1.4325 

1.3967 

Not .. 
Available 

0.2 

1.2999 

1.2819 

1.3667 

1.5315 

1.5046 

n 

1.4408 

1.4039 

t! 

O.U 

1-3095 

1.2911 

1.377U 

1-5371 

1-5059 

H 

1.4673 

1.4262 

Tf 

0.6 

1.3263 

1.3073 

1.3963 

1.5U73 

1.5078 

tt 

1.5164 

1.4662 

If 

0.8 

1.3517 

1.3326 

1.U253 

1.5647 

1.5099 

it 

1.6004 

1.5288 

It 

0.9 

1.3679 

1.3500 

1.UUU0 

1.5790 

1.5107 


1 . 6668 

1.5707 

TT 

1.2 

1.52U7 

- 

1.6330 

I.050U 

- 

tt 

2.0902 

- 

If 

l.U 

1.6lU0 

- 

1.7129 

0.8U07 

- 

tt 

2.0483 

- 

IT 

1.6 

1.6113 

- 

1.6973 

0. 6668 

- 

0.6885 

1.9493 

- 

2.0749 

1.8 

1.5732 

- 

1.6495 

O.508U 

- 

0.5055 

1.8400 

- 

1.9263 

2.0 

1.5225 

- 

1.5892 

0.3606 

- 

0.3244 

1.7200 

- 

1.8036 

2.2 

1.4583 

- 

1.5166 

0.02U2 

- 

0.0337 

1.5514 

- 

1 . 6989 

2 . U 

1.3859 

- 

l.UUoU 

0. 

- 

0. 

1.3515 

- 

1.6079 
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Table 4. Comparison of Predicted Lii 


Geometry 

M 

K 

P 

A = 0°, A = 3.5 

0.6 

3-792 

X = 1.0 


3-778 

A = 1+5°, A = 1.0 

0 . 

1.433 

A = 1.0 


1.431 

A = 45°, A = 2.0 

0 . 

2.263 

X = 1.0 


2.279 

A = 70°, A = 2.24 

2.0 

2.1173 

A = 0. 


Not available 

(Arrow Wing) 



A = 60°, A = 2.0 

1.5 

2.6060 

X = 0.5 


Not available 

(Cropped Arrow Wing) 



Model I 

0.2 

1.8928 

(Figure 10) 


1.8047 

Model II 

0.2 

1.3492 

(Figure 10 ) 


1.3353 

Model III 

0.2 

1.4596 

(Figure 10) 


1.4506 

Model IV 

0.2 

'1 

(Figure 10) 



Model V 

0.2 

1.9574 

(Figure 10 ) 


1.8948 


Factors for Various Configurations 


K 0 
v,Ze 

K 

v,se 

Methods 

2.481 

1.279 

Present 

2.480 

1.250 

ref. 9 

1.102 

2.489 

Present 

1.101 

2.412 

ref. 9 

2.040 

1.874 

Present 

2.037 

1-925 

ref. 9 

2.3656 

0 . 

Present 

2.3733 

0 . 

Exact 
(ref. 12 ) 

1.3788 

2.l46l 

Present 

1.4037 

2.262 

Exact 
(ref. 12 ) 

3-1123 

0 . 

Present 

3.1315 

0 . 

ref. 9 

1.9634 

0.7854 

Present 

1.8748 

0.7821 

ref. 9 

2.2297 

0.8665 

Present 

2.2898 

0.8803 

ref. 9 


1.8546 

Present 


1.7706 

ref . 9 

3-5684 I 

0 . 

. Present 

3.9653 | 

0 . 

ref. 9 







































Table 5* Comparison of Predicted Center-of-Force Factors for a Cropped Delta Wing. 

A = 0.874, A = 63 °, A = 0.4. (Reference Point at the Leading Edge of Root 
C 


M 


-x / c 
p r 

“Xfi /c 
te r 



-x /c 
se r 


Present 

ref. 9 

Present 

ref. 9 

Present 

ref. 9 

0 . 

0.4228 

0.4353 

0.3512 

0.3576 

0.8128 

0.8150 

0.2 

0.4226 

0.4355 

0.3519 

0.3582 

0.8126 

0 . 8l48 

0.4 

0.4217 

0.4361 

0.3544 

0.3601 

0.8119 

0.8144 

0.6 

0.4197 

0.4378 

0.3589 

0.3631 

0 . 8106 

0 . 8137 

0.8 

0.4156 

0.4428 

0 . 3668 

O .3658 

0.8082 

0 . 8129 

0.9 

0 . 4111 

0.4501 

0.3730 

0.3631 

0.806l 

0 . 8129 

- - - 

- “ - - - 



Exact 
(ref. 12) 

- - - - 

- - Exact 

(ref. 12) 

1.8 

0.5295 

Not Determined 

0.3932 

o.4o 

0.8082 

0.8095 

2.0 

0.5401 

?r 

0.3967 

0.40 

0.8107 

0.8095 

















0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 


€ 

Figure 4.- Predicted Pressure Distributions for Cropped Arrow 
Wing of A - 2.0, A = 0.5 at M = 1.5 and a = 1.0 
Radian . 










Figure 8.- Distribution of Leading Edge and Wing-Tip 

Suction Coefficients for Cropped Delta Wing 
of A = 0.874 and X = 0.4 at M = 0. 







3h 





Model I Model II Model III 

( A = 1-6) (A = 0.9514) (A = 1.076) 



Model IV 
(A = 0.2682) 



Model V 
(A = 1.8754) 


Figure 10.- Geometry of Family of Double Deltas of Reference 16 
and Referred in Table 4. 
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